Assessment of voltage harmonics’ impact on the maximum load capacity of the power supply transformer for the LHCD system

To study the influence of the voltage total harmonic distortion (THDV) and its spectrum on the harmonic loss factor (FHL) and the maximum allowable load capacity ratio (Imax (pu)) for the transformer of the LHCD system, we use MATLAB software to model the LHCD system and study the winding loss and maximum allowed load capacity of the LHCD system transformer when the supply voltage source is sinusoidal and non-sinusoidal, respectively. The calculation is carried out by the method of IEEE standard C57.110 and the method of considering the skin effect. The calculation results of both methods clearly show that the THDV value of the supply voltage has a significant effect on the harmonic dissipation factor (FHL) and the maximum allowable load capacity ratio (Imax (pu)) of the transformer. And the calculation method considering the skin effect increases the maximum allowable load capacity ratio of the LHCD system transformer by about 2%. The research results have important reference value for the future retrofit design of LHCD system transformers.

studies are all based on IEEE Standard C57.110 and UL Standard 1561, using the harmonic loss factor (F HL ) and K factor to determine the maximum allowable load capacity ratio (the max permissible rms non-sinusoidal load current under rated conditions (I max (pu)) of the transformer 10,14 .However, the basic assumption of IEEE Standard C57.110 is that winding eddy current losses for conductors smaller than 3 mm gradually increase with the square of the harmonic frequency, in line with the maximum winding losses.However, this assumption is wrong for conductors of larger dimensions above the fifth harmonic due to the skin effect 15,16 .
In 17 , the influence of harmonics on winding loss and temperature distribution characteristics was studied, but possible solutions or measures were ignored.In 18 , it shows that the load loss of transformers increases with higher harmonic distortion rates and frequencies, emphasizing the need for measures to mitigate the impact of harmonics on transformer performance.But this article does not discuss the potential impact on transformer lifespan, reliability, or overall system stability.
In addition, some studies have found that the supply voltage harmonics also affect the current harmonic distortion of nonlinear loads [19][20][21][22] , so the voltage harmonics may indirectly affect the transformer winding losses and the maximum allowable load capacity of the transformer.
In this paper, considering the skin effect, the F HL and I max (pu) values of dry-type transformers supplying LHCD systems under sinusoidal and non-sinusoidal voltages are studied, and a revised load loss calculation model for transformer harmonic currents is proposed.The aim is to analyze the indirect effects of voltage harmonics on transformer winding losses and maximum allowable load capacity.

Description of HVPS system
The 4.6 GHz/6 MW LHCD system has a total of four identical HVPS systems.For the convenience of illustration, we only select one HVPS system.Its main circuit configuration is shown in Fig. 1.It can be seen from Fig. 1 that the HVPS system is a complex system with a power supply system voltage of 10 kV and a power distribution cabinet.Other main components are a soft start switch cabinet, multi-winding rectifier transformer, PSM power supply, and control unit.The three-phase 10 kV AC power grid controls the incoming lines of two 2000 kVA rectifier transformers through a 10 kV power distribution cabinet with soft-start switchgear.
Each rectifier transformer provides power for 32 sets of identical PSM modules, and the output voltage of each PSM power module is 810VDC.In order to obtain higher DC output voltage, 64 groups of power modules are connected in series.By controlling the IGBTs of 64 PSM module power supplies, the maximum output voltage of the HVPS system can reach − 50 kV 3,7 .The main parameters of transformers are shown in Table 1.

Derating calculation method IEEE standard C57.110 derating method
According to the definition of IEEE standard C57.110 3.1, the value of transformer loss consists of three parts, which are no-load loss, load loss, and total loss.No-load loss is equal to excitation loss, load loss is equal to impedance loss, and total loss is equal to the sum of no-load loss and load loss.The load loss consists of I 2 R loss and stray loss.The Stray loss can be calculated by subtracting resistance loss from measured impedance loss.The stray losses are further divided into winding eddy-current losses (P EC ) and stray losses in components other than the windings (P OSL ).Then the total load loss can be expressed as Eq. ( 1 where P LL represents the load loss, P EC represents the winding eddy-current loss, P OSL represents the other stray loss, P is the I 2 R loss portion of the load loss.P can be calculated by multiplying the winding resistance R by the square of the total rms value of the load individual harmonic currents (I h ): In IEEE standard C57.110, to convenient perform the calculations in practice, Eq. ( 1) applied to rated load conditions is rewritten on a per-unit basis as follows in Eq. ( 3): where P LL-R (pu) is the per-unit load loss under rated conditions, P EC-R (pu) is the per-unit winding eddy-current loss under rated conditions, P OSL-R (pu) is the per-unit other stray loss under rated conditions.
For any defined non-sinusoidal load current, the winding eddy-current loss P EC can be expressed as: where h is the harmonic order, h max is the highest significant harmonic number, I h is the rms current at harmonic h (amperes), I R is the rms fundamental current under rated frequency and rated load conditions.The equation for the rms current in per-unit form (base current is rated current) under non-sinusoidal load currents can be expressed as: where I(pu) is the per-unit rms load current, I h (pu) is the per -unit rms current at harmonic h.
Then the rms value of the nonsinusoidal load current I can be given by: To convenient the calculation of the power supply capability of a transformer under balanced non-sinusoidal load conditions, IEEE standard C57.110 defines the harmonic loss indicators F HL and F HL-STR .
Then the P LL can be recalculated by combining Eqs. ( 1), ( 3), ( 4), ( 6), (7), and (8), as shown in Eq. ( 9): (1) For dry-type transformers, generally do not consider temperature rise due to other stray losses, since the heat generated is dissipated by cooling air.Then Eq. ( 9) can be recalculated by Eq. ( 10): Therefore, the maximum permissible rms non-sinusoidal load current (I max (pu)) for dry-type transformers under rated conditions or the ratio of the transformer's maximum permissible current (I max ) to the rated current (I R ) is given by Eq. ( 11):

Considering skin effect
The main limitation of Eq. ( 4) is that it only applies to conductors under 3 mm.However, However, when the conductor is large and the frequency flowing through the conductor is high, the resistance of the transformer winding will increase due to the skin effect and proximity effect 22 , resulting in a relatively high error in the calculation result of Eq. (4).To make up for the limitation of Eq. ( 4), a function is adopted from 23,24 .In this paper, this function is used to improve the harmonic factor for transformers with larger winding conductors, which is shown in Eq. ( 12).
where F( h ) is equal to 1 sinh −sin cosh − cos , = T σ ( T is the field strength, and σ is the depth of penetration), σ R is the depth of penetration at rated power frequency, is the depth of penetration under harmonic condition.
The maximum permissible rms non-sinusoidal load current considering skin effect can then be expressed as follows:

Simulation
The purpose of this study is to analyze the influence of Tokamak (EAST) voltage harmonics on F HL and I max I R values of the transformers supplying the HVPS.To this end, we built a system simulation model in MATLAB software.The simulation model includes a programmable voltage source, forty-three sets of the same PSM power modules (C = 3 mF, R = 700 Ω, L = 1.1 mH), two identical dry-type transformers, of which ratings 2000 kVA and 10 kV (wye)/(16 * 0.6 (delta) + 16 * 0.6 (wye)), the secondary winding consists of copper strands of 3.65 mm × 11 mm.Each transformer supplies power to 32 sets of PSM power modules.Note that in each experiment, the HVPS system will choose to invest in different numbers of PSM module groups according to the different objectives of each experiment.Since the parameters of the 64 groups of PSM modules are the same, the low-voltage windings of the transformer supplying power to the PSM modules have only two connection modes: wye and delta.Therefore, for the convenience of analysis, we only analyze two of the PSM modules (one powered by the wye low voltage winding and the other by the delta low voltage winding).

Considering skin effect
In this paper, the PSM module of the delta connection type is simulated under two voltage scenarios.The first scenario is under sinusoidal voltage conditions.The second scenario is under non-sinusoidal voltage conditions (voltage distortion from 0 to 15%).
Scenario 1 The voltage source is a sinusoidal voltage, and the simulation results in Table 2.The harmonic loss coefficient of winding eddy current F HL calculated by IEEE STANDARD C57.110 is 5.0608, and the transformer capacity I max (pu) is about 80.45% of its sinusoidal load current capacity.When considering the skin effect, the F * HL is 4.5372, and the I * max (pu) is about 82.35% of its sinusoidal load current capability.Scenario 2 In the simulation, under the condition that the power supply is a non-sinusoidal voltage (the seventh harmonic voltage content is 1% to 15% respectively), the values of F HL , F * HL , I max (pu) and I * max (pu) are calculated and recorded, and the results are shown in Fig. 2.

Wye-connected
In this paper, the PSM module of the wye connection type is simulated under two voltage scenarios.The first scenario is under sinusoidal voltage conditions.The second scenario is under non-sinusoidal voltage conditions (voltage distortion from 0 to 15%).
Scenario 1 The voltage source is a sinusoidal voltage, and the simulation results are tabled in Table 3. From the simulation results, it can be concluded that under the condition of sinusoidal voltage, the harmonic loss coefficient of winding eddy current F HL calculated by IEEE STANDARD C57.110 is 5.0428, and the transformer (10) P LL pu = I 2 pu 1 + F HL P EC−R pu    Table 3.The simulation results (F HL , F * HL , I max (pu) and I * max (pu)) of the wye -connected.In the simulation, under the condition that the power supply is a non-sinusoidal voltage (the seventh harmonic voltage content is 1% to 15% respectively), the values of F HL , F * HL , I max (pu) and I * max (pu) are calculated and recorded, and the results are shown in Fig. 3.

Discussion
For scenario 1: • The results (the wye-connected and the delta-connected) reveal that the percentage error between the eddy current loss factors F HL and F * HL calculated by the two methods is about 9.5% under the condition of sinusoidal voltage.
• The results (the wye-connected and the delta-connected) suggest that the maximum allowable RMS non- sinusoidal load current (I max (pu) and I * max (pu)) under rated conditions calculated by the two methods differs by about 2%.That is, the calculation method considering the skin effect can increase the maximum allowable rms non-sinusoidal load current under rated conditions by 2%.
• The current practice recommended in the IEEE std.C57 10-2018 standard for evaluating transformers for LHCD systems is conservative.
For scenario 2: • The results (wye-connected and delta-connected) of the eddy current loss factors F HL and F * HL show that when the voltage distortion rate increases from 1 to 15%, the eddy current loss factors F HL and F * HL calculated by the two methods also increase accordingly and start to increase rapidly from the voltage distortion rate of 3%.
The results (wye-connected and delta-connected) of the maximum allowable RMS non-sinusoidal load current (I max (pu) and I * max (pu)) show that when the voltage distortion rate increases from 1 to 15%, the maximum allowable RMS non-sinusoidal load current (I max (pu) and I * max (pu)) calculated by the two methods also decrease accordingly and start to decreases rapidly from the voltage distortion rate of 3%.

Conclusion
In this paper, the calculation model of transformer loss and maximum allowable current of LHCD system under voltage harmonic environment is established, and the functional relationship between the calculation method of IEEE standard C57.110 and the calculation method considering skin effect is deduced.On this basis, a modified transformer harmonic current load loss calculation model is proposed.Finally, the simulation verification is carried out with Simulink/MATLAB software, and the results show that: www.nature.com/scientificreports/(a) When the voltage harmonic distortion rate (THDV%) is greater than 3%, the calculation results of both methods show that the harmonic loss factor value of the LHCD system transformer increases significantly and the maximum allowable current value decreases significantly.(b) Compare the calculation results of the IEEE Standard C57.110 method with the calculation results of the method considering the skin effect: The harmonic loss factor H value calculated by the former is about 5% higher than that of the latter two.The maximum allowable current F value calculated by the former is nearly 2% lower than that of the latter.The results mean that in the high-frequency environment, the transformer winding loss of the LHCD system is significantly affected by the winding skin effect.
The research in this paper has reference significance for the analysis of transformer losses under the background of voltage harmonics, and can also provide theoretical and data support for the future optimal design and safe operation of LHCD system transformers.

Figure 2 .
Figure 2. The delta-connected.(a) The influence of THDV% on F HL and F * HL .(b) The influence of THDV% on I max (pu) and I * max (pu).

Figure 3 .
Figure 3.The wye-connected.(a) The influence of voltage distortion rate on F HL and F * HL .(b) The influence of voltage distortion rate on I max (pu) and I * max (pu).

Table 1 .
The electrical parameters of the transformer.

Table 2 .
The simulation results (F HL , F * HL , I max (pu) and I * max (pu) for the delta-connected.